Resonant magnetic X-ray scattering spectra in SDW Cr 

- ab initio study - 
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Using ab-initio band structure calculation based on the local density approximation, Cr K-edge 
resonant X-ray magnetic scattering spectra are analyzed in the spin density wave (SDW) state of 
chromium. We perform band structure calculation, assuming an ideal bcc lattice structure with the 
lattice constant observed at the spin-flip temperature Tsf and a commensurate SDW state with the 
propagation vector close to the observed value at Tsf. Taking account of the spin-orbit interaction, 
we obtain the orbital moment on each Cr site induced in proportion to the local spin moment, 
which is quite small, at most a tenth of those in nickel or iron. In spite of the tiny 3d orbital 
moment, the orbital polarization is found to have large fluctuations as a function of energy. We 
obtain the scattering intensity at the Cr K-edge on the SDW magnetic Bragg spot, which shows 
resonant enhancement in good agreement with the experiment. The 3d orbital polarization is found 
to be highly correlated with the intensity of the resonant main peak, indicating that the 4p orbital 
polarization is mainly induced by the 3d orbital polarization through the p-d hybridization. 
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I. INTRODUCTION 

Resonant X-ray diffraction technique has attracted 
much attention in the field of studying spin and or- 
bital properties of the 3d or 4/ electrons in transition 
metal or rear-earth metal compounds. Both resonant 
and non-resonant diffraction techniques have been re- 
vealing the charge, spin, orbital and lattice characters 
in a wide variety of compounds. In particular magnetic 
diffraction techniques have been applied to various mag- 
netic materials for studying the spin and orbital magnetic 
properties piiSiiiSiSiL2i2iiSiiiii2iiii4iiSii& In this paper we 
focus on the magnetic diffraction in the spin density wave 
in chromium 4iii 

Chromium is a prototypical itinerant antiferromagnet 
forming a spin density wave (SDW) state below the Neel 
temperature. A lot of studies have been carried out 
for revealing the electronic and magnetic properties of 
chromium 4iii*iii2iia Hill et aU studied the SDW state 
and the charge density wave (CDW) states by means of 
the X-ray scattering at the Cr i^-edge. They estimated 
the lattice distortion by analyzing the nonresonant inten- 
sity, but failed to observe resonant enhancement at the 
Cr K-edge. On the other hand Mannix et ali^i have suc- 
cessfully observed clear resonant enhancement at the Cr 
K-edge and analyzed the spectra using an atomic model. 
Noting that the resonant enhancement of the magnetic 
scattering intensity is mainly caused by the orbital po- 
larization (OP) in the 4p bandjfi^it is interesting that a 
clear resonant peak is observed in spite of the very small 
orbital moment in the ground stateiii 

We analyze the resonant X-ray magnetic scattering 
spectra in the SDW chromium, using an ab initio band 



structure calculation based on the local density approx- 
imation (LDA). We cannot rely on the tight-binding 
model for the 4p band, because the 4p wave function 
is greatly extended in space so that in the crystal most 
weight of the 4p wave function is distributed on the neigh- 
boring site or interstitial region. We first carry out a con- 
ventional band calculation without taking account of the 
spin-orbit interaction (SOI). We obtain the ground state 
potential which gives the SDW state, then calculate the 
wave function in conduction band under the potential. 
This calculation reproduces the same results as those al- 
ready reportedj^S*^ We take account of the SOI in the 
final stage. We find that the orbital moment is induced 
in proportion to the local spin moment on each site, al- 
though the orbital moment is as small as a tenth of those 
in ferromagnetic Fe, Co or WiM^ It should be noted 
here that the 3d OP fluctuates considerably as a func- 
tion of energy in spite of the tiny orbital moment in the 
ground state. 

Using the wave function thus obtained within the LDA, 
we calculate the scattering intensity. The calculation re- 
produces well the photon energy dependence on the SDW 
Bragg spot and the Fano-type dip in the lower energy side 
of the resonant main peakiii The Fano effect originates 
from the interference between the amplitudes of the reso- 
nant and nonresonant scattering processes. The intensity 
of the resonant main peak arises from the OP in the 4p 
band. In order to clarify how the Ap OP is induced, we 
calculate the spectra with turning on and off the SOI se- 
lectively in the 3d and 4p states. We find that the 4p OP 
is induced by the 3d OP through the p-d hybridization, 
and that the 3d OP strongly correlates with the resonant 
peak intensity. 



This reminds us of the mechanism of the magnetic 
circular dichroism (MCD) spectra of X-ray absorption 
at the transition metal K-edge. We studied the MCD 
spectra in Fe, Co, Ni22iS^ and MnaGaCSl, and revealed 
that the intensity of the dichroism spectra correlates 
with the 3d OP through the hybridization between the 
4p and 3d states iS2i22i2iiSi Also, the sensitivity to elec- 
tronic states at neighboring sites is closely related to the 
mechanism of the RXS intensity on the orbital order- 
ing Bragg spot in transition metal compounds. Detailed 
studie a^^i^7i^ffi^i?ii?fliF based on band structure calcula- 
tions have revealed that the intensity is mainly controlled 
by the lattice distortion through the hybridization of the 
transition-metal 4p states and the oxygen 2p states at 
neighboring sites. 

This paper is organized as follows. In the next section, 
we briefly describe the band structure calculation and 
the ground state. In Sec. Ill we formulate the magnetic 
scattering amplitude. We discuss the photon energy de- 
pendence of the scattering intensity and its relation to 
the electronic states in Sec. IV. The last section is de- 
voted to the concluding remarks. 



II. BAND STRUCTURE CALCULATION 

Chromium forms a bcc lattice structure. As cool- 
ing its magnetic state turns into a transverse spin den- 
sity wave (TSDW) state at its Neel temperature Tn = 
311 K and into a longitudinal spin density wave (LSDW) 
state at its spin-flip temperature Tsf — 122 K. In 
the TSDW (LSDW) state, magnetic moments are per- 
pendicular (parallel) to the SDW propagation vector. 
The wavelength of the SDW is incommensurate with 
the lattice periodicity.-^ Charge density wave (CDW) 
and lattice strain wave (LSW) are accompanied by the 

opTTr17.18.32..S3 

In the present calculation, we assume an ideal bcc lat- 
tice structure with the lattice constant a — 5.45ao, where 
ao = O.529A is Bohr radius. The SDW wavelength is as- 
sumed to be AsDW — 20a, which is very close to the 
observed value at the spin-flip temperature Tsf- Fig- 
ure Q] schematically depicts the assumed magnetic struc- 
ture. For convenience, the coordinate system is cho- 
sen such that the z direction is parallel to the magnetic 
moment. We flrst carry out the band calculation self- 
consistently without the SOI, using the Korringa-Kohn- 
Rostoker (KKR) method with muffln-tin (MT) approxi- 
mation. The size of the MT sphere is chosen such that 
neighboring spheres are touching each other. The maxi- 
mum number of k points n^ = nkx x nky x Ukz is 30 x 30 x 2 
in the first Brillouin zone (FBZ) for the LSDW state and 
16 X 2 X 16 for the TSDW state. Then we calculate the 
eigenvalues and eigenfunctions with adding the SOI term 
^-^V {r) izSz to the muffin-tin potential V{r). Selfcon- 
sistent iteration is not carried out. The one electron ex- 
cited states are calculated up to the energy ep + 2.0 ryd, 
where ep is the Fermi level. We neglect the spin flip 



terms £+S- and £-S+ in the SOI. We expect that the ne- 
glected terms cause only minor correction to the OP in 
transition metals and their compounds i2^»2i^ Since we 
are dealing with a large unit cell containing 40 inequiva- 
lent Cr atoms, this may help us to reduce the calculation 
volume. 
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FIG. 1: Schematics of SDW states with Qsdw = ^M- The 
arrows represent the direction and the magnitude of magnetic 
moment, (a) a simple AF states, (b) a TSDW state, a is the 
lattice parameter of simple bcc structure, (c) an LSDW state. 



We define the spin moment at the jth site as the dif- 
ference between electron numbers with spin up and spin 
down inside the MT sphere placed at the jth site. That 
carried by the valence d electrons /xf j^, is defined as 
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where Uj^ is the number of spin a electrons in the d 
state inside the jth MT sphere, and those carried by the 
valence s and p electrons are similarly defined. Another 
important quantity is the local orbital moment at each 
site. That carried by the valence d electrons at the jth 
site /ioj-bi i^ defined as 
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where m represents the magnetic quantum number, n'j™ 
is the number of the spin a electrons in the d state speci- 
fied by m inside the jth MT sphere. The quantity carried 
by the valence p electrons is similarly defined. 

Figure El shows the distribution of the local spin and 
orbital moments in the LSDW state. The contributions 
of the valence s and p electrons to the local spin mo- 
ments are negligible compared with those of the d elec- 
trons. The local orbital moments are also sinusoidally 
modulated, being proportional to the local spin moment. 
The magnitude of the local orbital moment is quite small. 
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FIG. 2: Spin(upper panel) and orbital(lower panel) moment 
distribution in LSDW state. The number of k points in the 
FBZ is 16 X 16 X 2. 
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FIG. 3: Partial density of states projected on the MT sphere 
at tq. Gray line and black line represent the DOS projected 
on the d and p states, respectively. The number of k points 
in the FBZ is 30 x 30 x 2. Fermi level ep is at OeV 



The maximum value is about 0.004/1, being smaller than 
a tenth of those in Fe, Co, and NiiSSiS^ Therefore the con- 
tribution of the orbital moment to the local moment is 
negligible. Note that the spin and orbital moment distri- 
bution in the TSDW state are nearly the same as those 
in the LSDW state, except for a slight reduction of the p 
orbital moment. 

The moments thus obtained may be written as 



Pj = Ml cos (QsDW • T, 
+M5 cos (5QsDW 



M3 cos (3QsDW 



(3) 



in terms of the amplitude of the fundamental wave 
Ml and the amplitude of the odd-order harmonics M3, 
M^,- ■ •, where Qsdw is a propagation vector of the SDW, 
QsDW = Qtsdw = ^ (0, §, 0) for the TSDW state and 
Qsdw = Qlsdw = ^ (0, 0, i§) for the LSDW state. Tj 
denotes the position vector of the jth site, tj = jaz (y) 



for n = 0,2,4, 



and Tj = jaz (y) -I- ^a (x + y + z) in 
the LSDW (TSDW) state. For the 3d spin moments, 
the coefficients in Eq. IpJ are evaluated as Mf'^ = 
0.393ft, M|^ = -0.026fi, Ml'^ = 0.0025ft in the LSDW 
state. These values are consistent with the reported band 
calculations^SiSl and the experiment pii 

The charge density is also modulated, but its magni- 
tude is underestimated compared with the experimental 
value4i If the LSW is taken into account, this drawback 
may be improved. 

Finally in this section, we discuss the partial density of 
states (PDOS) defined as the density of states projected 
on the p{£ = 1) or d{£ — 2) state inside the MT sphere 
placed at each site. Figure 13 shows those quantities at 
site To, which has the largest magnitude of spin moment. 
The PDOS in the TSDW state is nearly the same as that 
in the LSDW state. The PDOS in the previous band 



calculations! is limited in the energy range —7 ~ 5eV. 
The present result is essentially the same as them. The 
PDOS of the p state is significantly small compared with 
that of the d state. This is consistent with the fact that 
the wave function of the p state is greatly extended in 
space. 



III. MAGNETIC X-RAY SCATTERING 
SPECTRA 

A. Formulation 

1. nonresonant magnetic scattering 

We use the formula given by Blume and Gibbs^i^^ for 
evaluating the nonresonant magnetic scattering ampli- 
tude. The smallness of the orbital angular moment in 
Cr may allow to ignore the orbital contribution to the 
scattering amplitude. Then we obtain the expression of 
the amplitude /nr as 



/NR(q',e'^q,e) 
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where 



/j[^(q',e'^q,e)=A(-2* 
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Here, /^^ is the scattering amplitude of the atom at 
and TV is the number of unit cells in the crystal. Vectors q 
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and e represent the wave number and polarization of the 
incident photon, respectively, while q' and e' are those 
of the emitted photon. Unit vectors q and q' are defined 
as q/g and q'/g', respectively. Q is a scattering vector 
defined by q' — q. huj, ii^Ryd, and mc^ represent the 
photon energy, the energy constant equal to 1 Ryd, and 
the rest energy of an electron, respectively. Function 
G (Q) is defined such that it is unity when Q coincides 
with one of the reciprocal vectors for the SDW supper 
cell, and zero other than the case. Sz ' (Q) is a Fourier 
transform of the local spin moment inside the MT sphere 
at Tj . The direction of magnetization is assumed to be 
parallel to the z direction. The factor A is given by 



Within the spherical approximatioDSiS,, Sj (Q) is ex- 
pressed as 



Sl^ (Q) = 2 
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^{pT,Ar)-pT,i{r))dr, (10) 



where Rmt is the radius of the MT sphere and PTjc (r) 
is the distribution function of the electrons with spin a 
inside the MT sphere at Tj . With the help of the band 
calculation, we evaluate Eq. Ijl0|l and thereby Eqs. Q 
and ijnj. 
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where a 



he 



represents the fine structure constant. 




FIG. 4: Scattering geometry for the LSDW(left) and 
TSDW(right) configuration. Gray arrows represent the di- 
rection of the magnetic moment. 

In the following, we consider the scattering in the an' 
channel with the scattering vector Q = Qlsdw(tsdw) 
for the LSDW(TSDW) states. The scattering geometry 
is defined in Fig. 31 The scattering angle 9 and azimuthal 
angle ■0 are shown. The scattering angle is 6* « 21 deg. 
in both of the LSDW and TSDW states. We can eas- 
ily see that the scattering amplitude for the aa' channel 
vanishes in the LSDW state. The amplitude for the an' 
channel in the LSDW state is given by 



/nr (q', 



q, e) = A - 
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where 
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Here we ignore the photon energy dependence assuming 
huj « 5990 eV « 440ii;Ryd. Also the amplitude for an' 
channel in the TSDW state is given by 

/nr ~ -A^sin^ecosesmijSziCl). (9) 



2. resonant magnetic scattering 

We evaluate the resonant scattering amplitude exploit- 
ing the second order Fermi's golden rulejSiSi We have 
three terms, the orbital-orbital scattering term foo, the 
orbital-spin scattering term /os, and the spin-spin scat- 
tering term /ssi*>3-^ Since /os ~ ;^/oo ~ /oo x 10"^ 

and /ss ~ (;^)' /oo ~ /oo x 10"^ for the Cr if-edge, 
/os and /ss can be ignored in the present calculation. 
We consider only the dipole process, since the quadrupole 
process gives rise to only minor contribution. 

Accordingly the resonant scattering amplitude may be 
expressed as 



/R(q',e'^q,e) 

= 7VG(Q)^e^Q-^Vo6(q''e'^q,e), (11) 



where 
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and 
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Here, function 6 (x) is defined that it is unity when x > 
and zero other than the case. Yem (e) is spherical har- 
monics Yim{0,4>) with ex = sin^cosc?!), ey = sin sin 0, 
and Cz = cos 0. Wave function V'nka (r) is for the state 
specified by band index n, wave vector k, and spin a 
with energy e„ka; vector r represents the relative posi- 
tion from Tj. £isa represent the Is core level with spin a, 
and F is the core hole lifetime broadening. Wave function 



''pnka;Tjim (r) is the projection of ipnUa (r) on the state 
specified by azimuthal quantum number £ and magnetic 
quantum number m inside the MT sphere at tj (it is 
normalized inside the MT sphere). We write the wave 
function as 

i>nka;Tjim (r) = Rnka;T,ima (r) Yi„i (f ) , (14) 

where Rnka-Tjima (r) is the normalized radial wave func- 
tion. 

The dipole transition matrix element TnUa-Tjima is 
given by 



/•Rmt rf. 

= / Rn\^a-T jlma ijr) Rt Asa {r)r'^ dv, (15) 

Jo '^0 

where Rvjisa (r) is the radial wave function of the Is 
state with spin a at tj . The factor A' in Eq. lfT2|l is 



2 / \ 2 



E^yd 



(16) 



Since hcu « e„kQ, — ei^, this factor nearly equal to the 
factor A in the nonresonant scattering amplitude. Here- 
after, we replace A' by A. 

'J- . 

We rewrite the site scattering factor /qq as 
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where the dipole transition density matrix p,J„j (e) is 
defined by 
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Incidentally, the absorption intensity /abs in the dipole 
process may be given by 
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where 5A stands for the imaginary part of A. 



The matrix p^^„ (e) bears the symmetry compatible 



with the LSDW or TSDW states. Let p'^^ (e) be 
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in the LSDW(TSDW) state with Qlsdw(tsdw) = ^m- 
Note that the off-diagonal elements a, (3 and 7 are zero 
in the LSDW state, because of the four fold rotational 
symmetry along the z direction. In the TSDW state, this 
symmetry is not hold in a strict sense. However, one may 
expect that it is approximately hold, because the wave 
length of the SDW is quite long Asdw = 20a, leading to 
very small off-diagonal elements in the transition density 
matrix. 

The scattering amplitude for the resonant process may 
be rewritten as 

/R(q',e' ^q,e) 

^''-AiVG(Q)^rw(e')ri;,(e) 
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Then p (e, Q) has the form 
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From the above-mentioned symmetry argument on 
p'''^ (e), the off-diagonal elements 77, C are zero in the 
LSDW state. The scattering amplitude for the cttt' chan- 
nel in the LSDW state becomes 

/r (q', e"" ^ q, e'") = -i A7VG (Q) i sin 6 



de 



H^M) 



{fku — e — eisa + ir) /mc^ ' 



(25) 



It is easily seen that the amplitude for the aa' channel is 
zero. Also, in the TSDW state, the off-diagonal elements 
7], C are expected to be very small. Neglecting them, we 
obtain the scattering amplitude for the an' channel in 
the TSDW state as 
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The amplitude in the aa' is zero if the ofT-diagonal ele- 
ments rj, C are zero. 

Using the non-resonant and resonant scattering ampli- 
tude thus obtained, the scattering cross section can be 
written as, 



da 

In 



^^rl |/nr + /r|^ 
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where Ve 



represents the classical electron radius. 



IV. RESULTS AND DISCUSSION 

Absorption spectra 

At present, it is difficult to estimate accurately the Is 
core level eis and the core-hole lifetime broadening F from 
the calculation. We adjust these values such that the 
calculated K-eAge absorption intensity reproduces the 
photon energy dependence observed in the experiment iii 
We calculate the spectra using Eq. ltT9| . assuming that 
ei, = ep - 5988.4 eV and T = 0.6 - leV. Figure 
shows the calculated absorption spectra in comparison 
with the experiment. The calculated spectral shape is 
roughly proportional to the PDOS for the unoccupied 4p 
state. The calculated curves agree well with the experi- 
ment except for a shape around hw = 6010 eV. The peak 
around 6010 eV comes from the peak around 21 eV in 
the p PDOS. Note that no noticeable difference is found 
between the LSDW and TSDW states. This good agree- 
ment between the calculation and the experiment may 
justify the neglect of the core-hole potential in the calcu- 
lation of resonant scattering spectra. 
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FIG. 5: Comparison of the absorption spectra. Solid and 
broken lines represent the calculated spectra and experiment. 
Is core level is assumed to be eis = ep — 5988.4 eV. The 
lifetime broadening of the intermediate states is assumed to 
be r = 0.6eV and 1.0 eV. 



magnetic scattering spectra 

Putting F — 0.6 eV and using the same values of ei^ 
as in the absorption spectra, we calculate the magnetic 
scattering spectra. We first check the convergence of 
the spectral curve with varying number of fc-points rifc 
in the FBZ. Figure El shows the spectral curves calcu- 
lated with three values of nfc at Q = Qlsdw in the 
LSDW state. The spectra calculated at ftw > 6000 eV 
show strong dependence on n^, no sign of convergence 
even for n^ = 30 x 30 x 2. On the other hand, the curve 
around huj = 5990 eV depends little on nfc, suggesting a 
convergence even for nfc = 12 x 12 x 2. 
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FIG. 6: Convergence of the scattering intensity at Q = 
Qlsdw for LSDW state against the number of fc-points. 
The number of A;-points are 12 x 12 x 2, 16 x 16 x 2 and 
30 X 30 X 2 for LSDW(l), LSDW(2) and LSDW(3), respec- 
tively. The Is core level and the lifetime broadening assumed 
to be eis = ep — 5988.4 eV and F = 0.6 eV, respectively. 



Figure shows the calculated spectral curves in 
the LSDW and TSDW states, in comparison with the 
experiment^. The absorption correction is not made in 
the experimental curves. The calculation reproduces suc- 
cessfully not only the Fano-type dip around tujj = 5986 eV 
in the TSDW state but also the resonant behavior in good 
agreement with the experiment. 

The resonant intensity arises from the OP in the 4p 
bands. In order to clarify how the OP is induced, we 
investigate the change of the electronic states by turning 
on and off the SOI on the p and d states selectively. Fig- 
ure IHl shows the density of OP ^ij" (e) as a function of 
energy in the LSDW state. This quantity is defined as 
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where p]^ (e) is the partial DOS projected on the state 
specified by azimuthal quantum number £ and magnetic 
quantum number m inside the MT sphere at site tq. 
(The site Tq has the largest magnetic moment.) When 
the SOI is working on both the p and d states (panel a), 
/i2 " (e) for the d states is quite large for < e < 2 eV but 
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FIG. 7: (a) Calculated magnetic scattering intensities for (ttt' 
channel on the SDW Bragg spots Qtsdw in the TSDW state. 
The azimuthal angle is assumed to be i/j = 90deg. (b) Those 
on the spot Qlsdw in the LSDW state. The number of k 
points in FBZ is n^ = 16 x 16 x 2. The azimuthal angle is 
assumed to he ip = ^ for the TSDW state. Is core level 
is assumed to be eis = ep — 5988.4 eV. Both results with 
F = 0.6 eV and F = 1.0 eV are shown for comparison. The 
spectra with F = l.OeV is normalized so that the peak hight 
at 5990 eV is equal with that of the spectra with F — 0.6 eV. 
The experiment curve is traced from Refi--. 



nearly zero for e > 3 eV. fi[° (e) for the p states is smaller 
in two order of magnitude than //j ° i^)- When the d SOI 
is turned off but the p SOI is still turned on (panel b), 
fi^" (e) becomes very small in the whole energy range. 
Corresponding to this change, p^" (e) also becomes very 
small for < e < 5eV, although it does not change for 
e > 5 eV from the values in the case when both of the 
p and d SOI's are turned on. When only the p SOI is 
turned off but the d SOI is turned on (panel c), pj" (^) 
is essentially unchanged in the whole energy range from 
the value in the case when both of the p and d SOI's are 
turned on. ^'[° (e) changes little for e < 5eV, but it is 
considerably reduced for e > 10 eV. 

These results clearly indicate that the p OP for < 
e < 5 eV is mainly induced by the d OP through the 
p-d hybridization. Since the d DOS is concentrated in 
e < 2eV, the hybridization with the p states at neigh- 
boring site may be effective only for e < 2 eV. The effect 
of the p SOI on the p OP seems to be reduced by this hy- 
bridization. Note that the p and d states do not hybridize 
each other at the same site because the SDW wavelength 
is so long that each site almost keeps the inversion sym- 
metry, particularly at Tq. The p OP for e > 10 eV is 
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FIG. 8: Orbital moment densities p.-^ ' (e) and /ij ' (e) on the 
site at Tq. (a) Both of the p SOI and d SOI are turned on. 
(b) The d SOI is turned off. (c) The p SOI is turned off. 
Solid line and gray line represent the orbital density for the 
p state p,i ' (e) and d symmetric states /^j ^ (^) : respectively. 
The orbital density for the p states is multiplied by 100. The 
number of fc-points in FBZ is 30 x 30 x 2. 



mainly induced by the p SOI. 

Since the resonant intensity arises from the p OP, the 
above finding that the p OP for e < 5 eV is induced by 
the d OP through the p-d hybridization may be confirmed 
by calculating directly the spectra with turning on and 
off the p SOI and d SOI. Figure El shows the calculated 
results. When only the d SOI is turned off, the resonant 
main peak at huj = 5990 eV disappears, while the reso- 
nant structures for higher energies change little. On the 
other hand, when only the p SOI is turned off, the main 
peak dose not change, while the resonant structure for 
higher energies almost disappears. This SOI dependence 
of the spectra corresponds well to the SOI dependence of 
the p OP. Consequently, it is conclude that the resonant 
main peak at huj = 5990 eV strongly correlates with the 
3d states. 

Finally we discuss the origin of the Fano-type dip. 
Since the Fano effect arises from the interference between 
the non-resonant and resonant terms, we look into /nr 
and /r, which are shown in Fig. ^1 /nr is purely imag- 
inary. Note that /r in the LSDW state is nearly propor- 
tional to that in the TSDW state with reverse sign. The 
reverse sign arises from the difference of the scattering 
geometry in the LSDW and TSDW states. As a result, 
in the TSDW state, the imaginary part of /r and /nr 
can be canceled at fto; « 5987 eV, with a small real part 
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FIG. 9: Scattering intensity on the SDW spot Qlsdw in 
LSDW state, (a) The SOI in p state is turned on and the 
SOI in d state is turned off. (b) The SOI in p state is turned 
off and the SOI in d state in turned on. The Is core level is 
assumed as eis = e-p — 5988.4 eV. The gray line is shown for 
comparison, which is calculated with both of the p SOI and 
d SOI being turned on. 
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FIG. 10: (a) Scattering amplitudes on the spot Qtsdw for the 
TSDW state, (b) Those on the spot Qlsdw for the LSDW 
state. Solid and broken lines represent the real and imaginary 
part of the amplitude due to the resonant scattering, respec- 
tively. Thin broken line represents the imaginary part of the 
amplitude due to the nonresonant scattering. The real part 
of the nonresonant contributions is zero, eia — ep — 5988.4 eV 
andr= 1.0 eV. 



of /r. This is the reason why the Fano dip is seen at 
hw ~ 5987 eV. On the other hand, one can easily see 
that such a cancellation is hard in the LSDW state. 



V. CONCLUDING REMARKS 

We have performed an ab initio band structure calcu- 
lation based on the LDA with a large unit cell. Assum- 
ing that the potential has a modulation with Qsdw = 
^^, we reproduced well the SDW state previously 
reportedi2Si21 Since the orbital momentum is essential to 
the resonant magnetic scattering, we have taken account 
of the SOI within the izSz approximation. The orbital 
magnetic moment is mainly induced on the 3d states in 
proportion to the spin magnetic moment at each Cr site. 

We have calculated the absorption spectra at the Cr 
iiT-edge. The dipole matrix elements are evaluated by 
using the wave functions obtained from the band struc- 
ture calculation. Having neglected the core-hole poten- 
tial, we have obtained the absorption spectra in good 
agreement with the experiment Ji This may justify ne- 
glecting the core-hole potential in the calculation of the 
resonant magnetic scattering spectra. 

We have calculated the resonant x-ray magnetic scat- 
tering spectra on the SDW magnetic Bragg spot in the 
TSDW and LSDW states. A resonant enhancement has 



been reproduced near the ii'-edge in good agreement with 
the experiment »iii In the TSDW state, the off-diagonal 
matrix elements of the scattering amplitude have been 
neglected, since they are expected to be small from a sym- 
metry argument. Experiments on the azimuthal angle 
dependence of the resonant intensity may help us check 
the accuracy of the argument. 

The mechanism of the resonant magnetic scattering 
has been clarified. We found that the OP on the 3d states 
fluctuates considerably as a function of energy, although 
the the local orbital magnetic moment in the d states is 
quite small. This d OP induces the p OP around the 
Fermi level through the p-d hybridization, and thereby 
gives rise to the resonant scattering intensity. The res- 
onant main peak aX huj = 5990 eV strongly correlates 
with the 3d states. This finding is consistent with the 
analysis^! of the resonant magnetic scattering of UGa^i^, 
where a large signal is found at the Ga K-edge. This is 
interpreted as a refiection of the large OP on the uranium 
5/ states. 

We have reproduced a Fano-type dip at huj fn 5987 eV 
in the TSDW stateiii The Fano-type dip is seen only in 
the TSDW state. This difference arises from the sign 
of the resonant scattering amplitude (the energy depen- 
dence of the amplitude is nearly the same in the two 
states except for the sign) , owing to the different scatter- 



ing geometry. 

Finally we comment on the CDW state. The assump- 
tion of an ideal bcc lattice gives rise to only a small CDW 
amplitude. The LSW may be dispensable for obtaining a 
sufficient magnitude of the CDW amplitude. A detailed 
study on the CDW state is left in future. 
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